Circadian regulation of sinoatrial nodal cell pacemaking function: Dissecting the roles of autonomic control, body temperature, and local circadian rhythmicity.

Strong circadian (~24h) rhythms in heart rate (HR) are critical for flexible regulation of cardiac pacemaking function throughout the day. While this circadian flexibility in HR is sustained in diverse conditions, it declines with age, accompanied by reduced maximal HR performance. The intricate regulation of circadian HR involves the orchestration of the autonomic nervous system (ANS), circadian rhythms of body temperature (CRBT), and local circadian rhythmicity (LCR), which has not been fully understood. Here, we developed a mathematical model describing ANS, CRBT, and LCR in sinoatrial nodal cells (SANC) that accurately captures distinct circadian patterns in adult and aged mice. Our model underscores how the alliance among ANS, CRBT, and LCR achieves circadian flexibility to cover a wide range of firing rates in SANC, performance to achieve maximal firing rates, while preserving robustness to generate rhythmic firing patterns irrespective of external conditions. Specifically, while ANS dominates in promoting SANC flexibility and performance, CRBT and LCR act as primary and secondary boosters, respectively, to further enhance SANC flexibility and performance. Disruption of this alliance with age results in impaired SANC flexibility and performance, but not robustness. This unexpected outcome is primarily attributed to the age-related reduction in parasympathetic activities, which maintains SANC robustness while compromising flexibility. Our work sheds light on the critical alliance of ANS, CRBT, and LCR in regulating time-of-day cardiac pacemaking function and dysfunction, offering insights into novel therapeutic targets for the prevention and treatment of cardiac arrhythmias.

Mechanistic insight into the functional role of human sinoatrial node conduction pathways and pacemaker compartments heterogeneity: A computer model analysis.

The sinoatrial node (SAN), the primary pacemaker of the heart, is responsible for the initiation and robust regulation of sinus rhythm. 3D mapping studies of the ex-vivo human heart suggested that the robust regulation of sinus rhythm relies on specialized fibrotically-insulated pacemaker compartments (head, center and tail) with heterogeneous expressions of key ion channels and receptors. They also revealed up to five sinoatrial conduction pathways (SACPs), which electrically connect the SAN with neighboring right atrium (RA). To elucidate the role of these structural-molecular factors in the functional robustness of human SAN, we developed comprehensive biophysical computer models of the SAN based on 3D structural, functional and molecular mapping of ex-vivo human hearts. Our key finding is that the electrical insulation of the SAN except SACPs, the heterogeneous expression of If, INa currents and adenosine A1 receptors (A1R) across SAN pacemaker-conduction compartments are required to experimentally reproduce observed SAN activation patterns and important phenomena such as shifts of the leading pacemaker and preferential SACP. In particular, we found that the insulating border between the SAN and RA, is required for robust SAN function and protection from SAN arrest during adenosine challenge. The heterogeneity in the expression of A1R within the human SAN compartments underlies the direction of pacemaker shift and preferential SACPs in the presence of adenosine. Alterations of INa current and fibrotic remodelling in SACPs can significantly modulate SAN conduction and shift the preferential SACP/exit from SAN. Finally, we show that disease-induced fibrotic remodeling, INa suppression or increased adenosine make the human SAN vulnerable to pacing-induced exit blocks and reentrant arrhythmia. In summary, our computer model recapitulates the structural and functional features of the human SAN and can be a valuable tool for investigating mechanisms of SAN automaticity and conduction as well as SAN arrhythmia mechanisms under different pathophysiological conditions.

Spatially resolved multiomics of human cardiac niches.

The function of a cell is defined by its intrinsic characteristics and its niche: the tissue microenvironment in which it dwells. Here we combine single-cell and spatial transcriptomics data to discover cellular niches within eight regions of the human heart. We map cells to microanatomical locations and integrate knowledge-based and unsupervised structural annotations. We also profile the cells of the human cardiac conduction system1. The results revealed their distinctive repertoire of ion channels, G-protein-coupled receptors (GPCRs) and regulatory networks, and implicated FOXP2 in the pacemaker phenotype. We show that the sinoatrial node is compartmentalized, with a core of pacemaker cells, fibroblasts and glial cells supporting glutamatergic signalling. Using a custom CellPhoneDB.org module, we identify trans-synaptic pacemaker cell interactions with glia. We introduce a druggable target prediction tool, drug2cell, which leverages single-cell profiles and drug-target interactions to provide mechanistic insights into the chronotropic effects of drugs, including GLP-1 analogues. In the epicardium, we show enrichment of both IgG+ and IgA+ plasma cells forming immune niches that may contribute to infection defence. Overall, we provide new clarity to cardiac electro-anatomy and immunology, and our suite of computational approaches can be applied to other tissues and organs.

Lipopolysaccharide-induced sepsis impairs M2R-GIRK signaling in the mouse sinoatrial node.

Sepsis has emerged as a global health burden associated with multiple organ dysfunction and 20% mortality rate in patients. Numerous clinical studies over the past two decades have correlated the disease severity and mortality in septic patients with impaired heart rate variability (HRV), as a consequence of impaired chronotropic response of sinoatrial node (SAN) pacemaker activity to vagal/parasympathetic stimulation. However, the molecular mechanism(s) downstream to parasympathetic inputs have not been investigated yet in sepsis, particularly in the SAN. Based on electrocardiography, fluorescence Ca2+ imaging, electrophysiology, and protein assays from organ to subcellular level, we report that impaired muscarinic receptor subtype 2-G protein-activated inwardly-rectifying potassium channel (M2R-GIRK) signaling in a lipopolysaccharide-induced proxy septic mouse model plays a critical role in SAN pacemaking and HRV. The parasympathetic responses to a muscarinic agonist, namely IKACh activation in SAN cells, reduction in Ca2+ mobilization of SAN tissues, lowering of heart rate and increase in HRV, were profoundly attenuated upon lipopolysaccharide-induced sepsis. These functional alterations manifested as a direct consequence of reduced expression of key ion-channel components (GIRK1, GIRK4, and M2R) in the mouse SAN tissues and cells, which was further evident in the human right atrial appendages of septic patients and likely not mediated by the common proinflammatory cytokines elevated in sepsis.

Non-canonical role of the sympathetic nervous system in the day-night rhythm in heart rate.

Although, for many decades, the day-night rhythm in resting heart rate has been attributed to the parasympathetic branch of the autonomic nervous system (high vagal tone during sleep), recently we have shown that there is a circadian clock in the cardiac pacemaker, the sinus node, and the day-night rhythm in heart rate involves an intrinsic rhythmic transcriptional remodelling of pacemaker ion channels, particularly Hcn4. We have now investigated the role of the sympathetic branch of the autonomic nervous system in this and shown it to have a non-canonical role. In mice, sustained long-term block of cardiac ß-adrenergic receptors by propranolol administered in the drinking water abolished the day-night rhythm in pacemaking in the isolated sinus node. Concomitant with this, there was a loss of the normal day-night rhythm in many pacemaker ion channel transcripts. However, there was little or no change in the local circadian clock, indicating that the well-known day-night rhythm in sympathetic nerve activity is directly involved in pacemaker ion channel transcription. The day-night rhythm in pacemaking helps explain the occurrence of clinically significant bradyarrhythmias during sleep, and this study improves our understanding of this pathology. This article is part of the theme issue 'The heartbeat: its molecular basis and physiological mechanisms'.

Genetics of sinoatrial node function and heart rate disorders.

The sinoatrial node (SAN) is the primary pacemaker of the mammalian heart, initiating its electrical activation and ensuring that the heart's functional cardiac output meets physiological demand. SAN dysfunction (SND) can cause complex cardiac arrhythmias that can manifest as severe sinus bradycardia, sinus arrest, chronotropic incompetence and increased susceptibility to atrial fibrillation, among other cardiac conditions. SND has a complex aetiology, with both pre-existing disease and heritable genetic variation predisposing individuals to this pathology. In this Review, we summarize the current understanding of the genetic contributions to SND and the insights that they provide into this disorder's underlying molecular mechanisms. With an improved understanding of these molecular mechanisms, we can improve treatment options for SND patients and develop new therapeutics.

cAMP signaling affects age-associated deterioration of pacemaker beating interval dynamics.

Sinoatrial node (SAN) beating interval variability (BIV) and the average beating interval (BI) are regulated by a coupled-clock system, driven by Ca2+-calmodulin activated adenylyl cyclase, cAMP, and downstream PKA signaling. Reduced responsiveness of the BI and BIV to submaximal, [X]50, ß-adrenergic receptor (ß-AR) stimulation, and phosphodiesterase inhibition (PDEI) have been documented in aged SAN tissue, whereas the maximal responses, [X]max, do not differ by age. To determine whether age-associated dysfunction in cAMP signaling leads to altered responsiveness of BI and BIV, we measured cAMP levels and BI in adult (2-4 months n = 27) and aged (22-26 months n = 25) C57/BL6 mouse SAN tissue in control and in response to ß-AR or PDEI at X50 and [X]max. Both cAMP and average BI in adult SAN were reduced at X50, whereas cAMP and BI at Xmax did not differ by age. cAMP levels and average BI were correlated both within and between adult and aged SAN. BIV parameters in long- and short-range terms were correlated with cAMP levels for adult SAN. However, due to reduced cAMP within aged tissues at [X]50, these correlations were diminished in advanced age. Thus, cAMP level generated by the coupled clock mechanisms is tightly linked to average BI. Reduced cAMP level at X50 in aged SAN explains the reduced responsiveness of the BI and BIV to ß-AR stimulation and PDEI.

Local tissue mechanics control cardiac pacemaker cell embryonic patterning.

Cardiac pacemaker cells (CPCs) initiate the electric impulses that drive the rhythmic beating of the heart. CPCs reside in a heterogeneous, ECM-rich microenvironment termed the sinoatrial node (SAN). Surprisingly, little is known regarding the biochemical composition or mechanical properties of the SAN, and how the unique structural characteristics present in this region of the heart influence CPC function remains poorly understood. Here, we have identified that SAN development involves the construction of a "soft" macromolecular ECM that specifically encapsulates CPCs. In addition, we demonstrate that subjecting embryonic CPCs to substrate stiffnesses higher than those measured in vivo results in loss of coherent electrical oscillation and dysregulation of the HCN4 and NCX1 ion channels required for CPC automaticity. Collectively, these data indicate that local mechanics play a critical role in maintaining the embryonic CPC function while also quantitatively defining the range of material properties that are optimal for embryonic CPC maturation.

Development of a new histological identification method of human sinoatrial node suitable for immunohistochemical study.

Histological identification of the human sinoatrial node (SAN) remains a challenge. Conventional identification methods, such as Lev's method, have certain limitations. The aim of our study was to develop a new histological identification method that could properly identify the sinoatrial node, applicable to the immunohistochemical study of intra-nodal structures. Thirty-nine human autopsied hearts were included in this study. The cases included 23 men and 16 women ranging in age from 20 to 99 years. The sinoatrial area from eight control samples was cut in the vertical section using the conventional Lev's method. In our new method, called the "En face one-block method," the sinoatrial node was cut in "En face" at the junction of the right border of the right appendage and superior vena cava, placed in one long cassette, and serially cut using a microtome. Immunostaining was performed using primary antibodies against CD31, podoplanin (D2-40), S-100, and other proteins. The average area of the SAN on the slide glass in our new method was 32.2 mm2, which was significantly larger than that (3.59 mm2) of the control samples by Lev's method. The SAN area was positively correlated with age (r = 0.357; p = 0.026), especially in women (r = 0.626; p = 0.0095). The SAN group had significantly lower percentage of CD31-positive blood capillaries, higher percentage of podoplanin-positive lymphatic channels, and S-100-positive peripheral nerves. We successfully developed a novel cutting method applicable to immunohistochemical studies, with which we could provide a bird's-eye view of the sinoatrial nodes.

Pacemaker activity and ion channels in the sinoatrial node cells: MicroRNAs and arrhythmia.

The primary pacemaking activity of the heart is determined by a spontaneous action potential (AP) within sinoatrial node (SAN) cells. This unique AP generation relies on two mechanisms: membrane clocks and calcium clocks. Nonhomologous arrhythmias are caused by several functional and structural changes in the myocardium. MicroRNAs (miRNAs) are essential regulators of gene expression in cardiomyocytes. These miRNAs play a vital role in regulating the stability of cardiac conduction and in the remodeling process that leads to arrhythmias. Although it remains unclear how miRNAs regulate the expression and function of ion channels in the heart, these regulatory mechanisms may support the development of emerging therapies. This study discusses the spread and generation of AP in the SAN as well as the regulation of miRNAs and individual ion channels. Arrhythmogenicity studies on ion channels will provide a research basis for miRNA modulation as a new therapeutic target.

cAMP-PKA signaling modulates the automaticity of human iPSC-derived cardiomyocytes.

Human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) have been used to screen and characterize drugs and to reveal mechanisms underlying cardiac diseases. However, before hiPSC-CMs can be used as a reliable experimental model, the physiological mechanisms underlying their normal function should be further explored. Accordingly, a major feature of hiPSC-CMs is automaticity, which is regulated by both Ca2+ and membrane clocks. To investigate the mechanisms coupling these clocks, we tested three hypotheses: (1) normal automaticity of spontaneously beating hiPSC-CMs is regulated by local Ca2+ releases (LCRs) and cAMP/PKA-dependent coupling of Ca2+ clock to M clock; (2) the LCR period indicates the level of crosstalk within the coupled-clock system; and (3) perturbing the activity of even one clock can lead to hiPSC-CM-altered automaticity due to diminished crosstalk within the coupled-clock system. By measuring the local and global Ca2+ transients, we found that the LCRs properties are correlated with the spontaneous beat interval. Changes in cAMP-dependent coupling of the Ca2+ and M clocks, caused by a pharmacological intervention that either activates the ß-adrenergic or cholinergic receptor or upregulates/downregulates PKA signaling, affected LCR properties, which in turn altered hiPSC-CMs automaticity. Clocks' uncoupling by attenuating the pacemaker current If or the sarcoplasmic reticulum Ca2+ kinetics, decreased hiPSC-CMs beating rate, and prolonged the LCR period. Finally, LCR characteristics of spontaneously beating (at comparable rates) hiPSC-CMs and rabbit SAN are similar. In conclusion, hiPSC-CM automaticity is controlled by the coupled-clock system whose function is mediated by Ca2+-cAMP-PKA signaling.

The funny current If is essential for the fight-or-flight response in cardiac pacemaker cells.

The sympathetic nervous system fight-or-flight response is characterized by a rapid increase in heart rate, which is mediated by an increase in the spontaneous action potential (AP) firing rate of pacemaker cells in the sinoatrial node. Sympathetic neurons stimulate sinoatrial myocytes (SAMs) by activating ß adrenergic receptors (ßARs) and increasing cAMP. The funny current (If) is among the cAMP-sensitive currents in SAMs. If is critical for pacemaker activity, however, its role in the fight-or-flight response remains controversial. In this study, we used AP waveform analysis, machine learning, and dynamic clamp experiments in acutely isolated SAMs from mice to quantitatively define the AP waveform changes and role of If in the fight-or-flight increase in AP firing rate. We found that while ßAR stimulation significantly altered nearly all AP waveform parameters, the increase in firing rate was only correlated with changes in a subset of parameters (diastolic duration, late AP duration, and diastolic depolarization rate). Dynamic clamp injection of the ßAR-sensitive component of If showed that it accounts for ∼41% of the fight-or-flight increase in AP firing rate and 60% of the decrease in the interval between APs. Thus, If is an essential contributor to the fight-or-flight increase in heart rate.

Disorder in Ca2+ release unit locations confers robustness but cuts flexibility of heart pacemaking.

Excitation-contraction coupling kinetics is dictated by the action potential rate of sinoatrial-nodal cells. These cells generate local Ca releases (LCRs) that activate Na/Ca exchanger current, which accelerates diastolic depolarization and determines the pace. LCRs are generated by clusters of ryanodine receptors, Ca release units (CRUs), residing in the sarcoplasmic reticulum. While CRU distribution exhibits substantial heterogeneity, its functional importance remains unknown. Using numerical modeling, here we show that with a square lattice distribution of CRUs, Ca-induced-Ca-release propagation during diastolic depolarization is insufficient for pacemaking within a broad range of realistic ICaL densities. Allowing each CRU to deviate randomly from its lattice position allows sparks to propagate, as observed experimentally. As disorder increases, the CRU distribution exhibits larger empty spaces and simultaneously CRU clusters, as in Poisson clumping. Propagating within the clusters, Ca release becomes synchronized, increasing action potential rate and reviving pacemaker function of dormant/nonfiring cells. However, cells with fully disordered CRU positions could not reach low firing rates and their ß-adrenergic-receptor stimulation effect was substantially decreased. Inclusion of Cav1.3, a low-voltage activation L-type Ca channel isoform into ICaL, strongly increases recruitment of CRUs to fire during diastolic depolarization, increasing robustness of pacemaking and complementing effects of CRU distribution. Thus, order/disorder in CRU locations along with Cav1.3 expression regulates pacemaker function via synchronization of CRU firing. Excessive CRU disorder and/or overexpression of Cav1.3 boosts pacemaker function in the basal state, but limits the rate range, which may contribute to heart rate range decline with age and disease.

Structural and Electrical Remodeling of the Sinoatrial Node in Diabetes: New Dimensions and Perspectives.

The sinoatrial node (SAN) is composed of highly specialized cells that mandate the spontaneous beating of the heart through self-generation of an action potential (AP). Despite this automaticity, the SAN is under the modulation of the autonomic nervous system (ANS). In diabetes mellitus (DM), heart rate variability (HRV) manifests as a hallmark of diabetic cardiomyopathy. This is paralleled by an impaired regulation of the ANS, and by a pathological remodeling of the pacemaker structure and function. The direct effect of diabetes on the molecular signatures underscoring this pathology remains ill-defined. The recent focus on the electrical currents of the SAN in diabetes revealed a repressed firing rate of the AP and an elongation of its tracing, along with conduction abnormalities and contractile failure. These changes are blamed on the decreased expression of ion transporters and cell-cell communication ports at the SAN (i.e., HCN4, calcium and potassium channels, connexins 40, 45, and 46) which further promotes arrhythmias. Molecular analysis crystallized the RGS4 (regulator of potassium currents), mitochondrial thioredoxin-2 (reactive oxygen species; ROS scavenger), and the calcium-dependent calmodulin kinase II (CaMKII) as metabolic culprits of relaying the pathological remodeling of the SAN cells (SANCs) structure and function. A special attention is given to the oxidation of CaMKII and the generation of ROS that induce cell damage and apoptosis of diabetic SANCs. Consequently, the diabetic SAN contains a reduced number of cells with significant infiltration of fibrotic tissues that further delay the conduction of the AP between the SANCs. Failure of a genuine generation of AP and conduction of their derivative waves to the neighboring atrial myocardium may also occur as a result of the anti-diabetic regiment (both acute and/or chronic treatments). All together, these changes pose a challenge in the field of cardiology and call for further investigations to understand the etiology of the structural/functional remodeling of the SANCs in diabetes. Such an understanding may lead to more adequate therapies that can optimize glycemic control and improve health-related outcomes in patients with diabetes.

Emergence of heartbeat frailty in advanced age I: perspectives from life-long EKG recordings in adult mice.

The combined influences of sinoatrial nodal (SAN) pacemaker cell automaticity and its response to autonomic input determine the heart's beating interval variability and mean beating rate. To determine the intrinsic SAN and autonomic signatures buried within EKG RR interval time series change in advanced age, we measured RR interval variability before and during double autonomic blockade at 3-month intervals from 6 months of age until the end of life in long-lived (those that achieved the total cohort median life span of 24 months and beyond) C57/BL6 mice. Prior to 21 months of age, time-dependent changes in intrinsic RR interval variability and mean RR interval were relatively minor. Between 21 and 30 months of age, however, marked changes emerged in intrinsic SAN RR interval variability signatures, pointing to a reduction in the kinetics of pacemaker clock mechanisms, leading to reduced synchronization of molecular functions within and among SAN cells. This loss of high-frequency signal processing within intrinsic SAN signatures resulted in a marked increase in the mean intrinsic RR interval. The impact of autonomic signatures on RR interval variability were net sympathetic and partially compensated for the reduced kinetics of the intrinsic SAN RR interval variability signatures, and partially, but not completely, shifted the EKG RR time series intervals to a more youthful pattern. Cross-sectional analyses of other subsets of C57/BL6 ages indicated that at or beyond the median life span of our longitudinal cohort, noncardiac, constitutional, whole-body frailty was increased, energetic efficiency was reduced, and the respiratory exchange ratio increased. We interpret the progressive reduction in kinetics in intrinsic SAN RR interval variability signatures in this context of whole-body frailty beyond 21 months of age to be a manifestation of "heartbeat frailty."

Pacemaking function of two simplified cell models.

Simplified nonlinear models of biological cells are widely used in computational electrophysiology. The models reproduce qualitatively many of the characteristics of various organs, such as the heart, brain, and intestine. In contrast to complex cellular ion-channel models, the simplified models usually contain a small number of variables and parameters, which facilitates nonlinear analysis and reduces computational load. In this paper, we consider pacemaking variants of the Aliev-Panfilov and Corrado two-variable excitable cell models. We conducted a numerical simulation study of these models and investigated the main nonlinear dynamic features of both isolated cells and 1D coupled pacemaker-excitable systems. Simulations of the 2D sinoatrial node and 3D intestine tissue as application examples of combined pacemaker-excitable systems demonstrated results similar to obtained previously. The uniform formulation for the conventional excitable cell models and proposed pacemaker models allows a convenient and easy implementation for the construction of personalized physiological models, inverse tissue modeling, and development of real-time simulation systems for various organs that contain both pacemaker and excitable cells.

In vivo and ex vivo electrophysiological study of the mouse heart to characterize the cardiac conduction system, including atrial and ventricular vulnerability.

The mouse is a common and cost-effective animal model for basic research, and the number of genetically engineered mouse models with cardiac phenotype is increasing. In vivo electrophysiological study in mice is similar to that performed in humans. It is indispensable for acquiring intracardiac electrocardiogram recordings and determining baseline cardiac cycle intervals. Furthermore, the use of programmed electrical stimulation enables determination of parameters such as sinoatrial conduction time, sinus node recovery time, atrioventricular-nodal conduction properties, Wenckebach periodicity, refractory periods and arrhythmia vulnerability. This protocol describes specific procedures for determining these parameters that were adapted from analogous human protocols for use in mice. We include details of ex vivo electrophysiological study, which provides detailed insights into intrinsic cardiac electrophysiology without external influences from humoral and neural factors. In addition, we describe a heart preparation with intact innervation by the vagus nerve that can be used as an ex vivo model for vagal control of the cardiac conduction system. Data acquisition for in vivo and ex vivo electrophysiological study takes ~1 h per mouse, depending on the number of stimulation protocols applied during the procedure. The technique yields highly reliable results and can be used for phenotyping of cardiac disease models, elucidating disease mechanisms and confirming functional improvements in gene therapy approaches as well as for drug and toxicity testing.

Loss of Natriuretic Peptide Receptor C Enhances Sinoatrial Node Dysfunction in Aging and Frail Mice.

Heart rate (HR) is controlled by the sinoatrial node (SAN). SAN dysfunction is highly prevalent in aging; however, not all individuals age at the same rate. Rather, health status during aging is affected by frailty. Natriuretic peptides regulate SAN function in part by activating natriuretic peptide receptor C (NPR-C). The impacts of NPR-C on HR and SAN function in aging and as a function of frailty are unknown. Frailty was measured in aging wild-type and NPR-C knockout (NPR-C-/-) mice using a mouse clinical frailty index (FI). HR and SAN structure and function were investigated using intracardiac electrophysiology in anesthetized mice, high-resolution optical mapping in intact atrial preparations, histology, and molecular biology. NPR-C-/- mice rapidly became frail leading to shortened life span. HR was reduced and SAN recovery time was increased in older versus younger mice, and these changes were exacerbated in NPR-C-/- mice; however, there was substantial variability among age groups and genotypes. HR and SAN recovery time were correlated with FI score and fell along a continuum regardless of age or genotype. Optical mapping demonstrates impairments in SAN function that were also correlated with FI score. SAN fibrosis was increased in aged and NPR-C-/- mice and was graded by FI score. Loss of NPR-C results in accelerated aging and rapid decline in health status in association with impairments in HR and SAN function. Frailty assessment was effective and better able to distinguish aging-dependent changes in SAN function in the setting of shortened life span due to loss of NPR-C.

Slowing down as we age: aging of the cardiac pacemaker's neural control.

The cardiac pacemaker ignites and coordinates the contraction of the whole heart, uninterruptedly, throughout our entire life. Pacemaker rate is constantly tuned by the autonomous nervous system to maintain body homeostasis. Sympathetic and parasympathetic terminals act over the pacemaker cells as the accelerator and the brake pedals, increasing or reducing the firing rate of pacemaker cells to match physiological demands. Despite the remarkable reliability of this tissue, the pacemaker is not exempt from the detrimental effects of aging. Mammals experience a natural and continuous decrease in the pacemaker rate throughout the entire lifespan. Why the pacemaker rhythm slows with age is poorly understood. Neural control of the pacemaker is remodeled from birth to adulthood, with strong evidence of age-related dysfunction that leads to a downshift of the pacemaker. Such evidence includes remodeling of pacemaker tissue architecture, alterations in the innervation, changes in the sympathetic acceleration and the parasympathetic deceleration, and alterations in the responsiveness of pacemaker cells to adrenergic and cholinergic modulation. In this review, we revisit the main evidence on the neural control of the pacemaker at the tissue and cellular level and the effects of aging on shaping this neural control.

ß-Adrenergic Stimulation Synchronizes a Broad Spectrum of Action Potential Firing Rates of Cardiac Pacemaker Cells toward a Higher Population Average.

The heartbeat is initiated by pacemaker cells residing in the sinoatrial node (SAN). SAN cells generate spontaneous action potentials (APs), i.e., normal automaticity. The sympathetic nervous system increases the heart rate commensurate with the cardiac output demand via stimulation of SAN ß-adrenergic receptors (ßAR). While SAN cells reportedly represent a highly heterogeneous cell population, the current dogma is that, in response to ßAR stimulation, all cells increase their spontaneous AP firing rate in a similar fashion. The aim of the present study was to investigate the cell-to-cell variability in the responses of a large population of SAN cells. We measured the ßAR responses among 166 single SAN cells isolated from 33 guinea pig hearts. In contrast to the current dogma, the SAN cell responses to ßAR stimulation substantially varied. In each cell, changes in the AP cycle length were highly correlated (R2 = 0.97) with the AP cycle length before ßAR stimulation. While, as expected, on average, the cells increased their pacemaker rate, greater responses were observed in cells with slower basal rates, and vice versa: cells with higher basal rates showed smaller responses, no responses, or even decreased their rate. Thus, ßAR stimulation synchronized the operation of the SAN cell population toward a higher average rate, rather than uniformly shifting the rate in each cell, creating a new paradigm of ßAR-driven fight-or-flight responses among individual pacemaker cells.